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ABSTRACT 

This is a report on detailed modeling of young high-mass protostellar candidates during their most 
embedded and obscured phases. We performed narrowband mid-infrared imaging of three candidate 
high-mass protostellar objects in Gil. 94-0. 62, G29.96-0.02, and G45. 07-^0. 13 at Gemini Observatory 
using the Thermal- Region Camera and Spectrograph (T-ReCS). The sources were imaged through up 
to 11 narrowband filters, sampling their SEDs over the entire 2-25 /im infrared range. For the first 
time, we have fit the observed SEDs of massive protostars with models that take into account departures 
from spherical symmetry in the infalling envelopes. In this way, we have been able to back out of the 
models detailed physical parameters for these earliest stages of massive stellar life. Our detailed modeling 
suggests that massive star formation can proceed in a way very similar to the formation of low-mass 
stars. 

Subject headings: circumstellar matter — infrared: ISM — stars: formation — ISM: individual 
(Gll.94-0.62, G29.96-0.02, G45.07-F0.13) — stars: early type 



1. INTRODUCTION 

The most massive stars begin their lives in the darkest 
recesses of giant molecular clouds, obscured from our view 
by a thick veil of dust and gas. Most of the time the ear- 
liest stages of massive star formation cannot be observed 
at visual wavelengths, and in many cases even the near- 
infrared emission from these protostars cannot escape. We 
owe much of our ignorance of the massive star formation 
process to the fact that we have a very limited range of 
wavelengths at our disposal with which to observe these 
youngest stages of massive stellar evolution. 

The problem is further compounded by the fact that 
massive OB stars (> 8 M©) are relatively few in number 
and, like all stages of massive stellar evolution, the forma- 
tion process proceeds very quickly. It is therefore difficult 
to catch massive stars in the act of forming. Addition- 
ally, there is the problem of observing these objects with 
sufficient angular resolution. Massive star-forming regions 
generally lie at kiloparsec distances, making them more 
than an order of magnitude farther away than the closest 
low-mass star-forming regions. Furthermore, the percent- 
age of binaries and multiple star systems increases as a 
function of mass (Larson 2001), and consequently it is be- 
lieved that the vast majority of massive stars do not form 
alone. It is this combined distance and multiplicity prob- 
lem that makes it difficult to resolve and study in detail a 
single massive star in its earliest years of formation. 

Observations that have been extremely helpful in the 
study of massive stars have been cm- wavelength radio con- 
tinuum and molecular line emission observations. Interfer- 



ometric arrays can observe with subarcsecond angular res- 
olution the radio continuum emission from the photoion- 
ized region in the near-stellar environment around young 
massive stars. These compact or ultra-compact H II (UC 
H II) regions are thought to exist ~10^ years after the 
formation of a massive star (De Free, Rodriguez, & Goss 
1995). Until recently, this was the earliest observed phase 
of massive star formation. 

In order to understand better their formation process, 
we need to observe massive stars at ages even younger 
than 10^ years. This means trying to observe sources be- 
fore they develop UC H II regions, which in most of the 
cases limits our window of continuum observations to a 
relatively narrow wavelength regime from the mid-infrared 
(~3 /im) to the millimeter (~1 mm). Observing sources 
in this "high mass protostellar object" (HMPO) phase 
presents a challenge since the best resolutions achievable 
throughout most of this wavelength range (namely the far- 
infrared to the sub-millimeter) have, until recently, been 
on the order of a few to tens of arcseconds. However, the 
bright molecular line emission from these sources would 
still be observable with high angular resolution radio in- 
terferometers. 

Hot cores of compact and dense molecular gas have been 
known to exist in massive star forming regions for sev- 
eral decades now (e.g., Wilson, Downes, & Bieging 1979). 
Some of these "hot molecular cores" (HMCs) are sim- 
ply externally heated knots of material, but others may 
have stars forming at their centers. Recently some HMCs 
have been gaining serious attention as possibly forming 
extremely young massive stars at their centers (Cesaroni 
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et al. 1994; Kurtz et al. 2000). Walmsley (1995) proposed 
that these types of HMCs are undergoing an intense accre- 
tion phase which prevents the development of a UC H II 
region. Since these types of HMCs do not possess radio 
continuum emission and arc internally heated by massive 
protostars, they fit into the category of HMPOs. 

Though the compositions of these sources are begin- 
ning to be understood in detail through the molecular line 
studies (e.g., Rodgers & Charnley 2001), continuum ob- 
servations remain few because of the angular resolution 
problem. However, continuum observations from the mid- 
infrared to millimeter have the potential for revealing valu- 
able information about the massive star formation process. 
Models of the spectral energy distribution (SED) of HMCs 
have been developed under the assumption that they have 
a central massive star undergoing spherical accretion of a 
free-falling envelope of dust and gas (Osorio, Lizano, & 
D'Alessio 1999). The success of these models in fitting 
the observable characteristics of HMCs gave strong sup- 
port to the hypothesis that some of these objects were 
true HMPOs (i.e., not starless HMCs) in a stage previous 
to the development of an UC H II region. However, for 
that work, only mid-infrared upper limits were available 
in most of the sources. Since the mid-infrared range of the 
SED is highly sensitive to the geometry of the source (i.e., 
degree of flattening), in principle, these kinds of observa- 
tions would allow testing of more sophisticated models, 
with a degree of detail similar to that of those that have 
been developed for low-mass protostars (e.g., Osorio et al. 
2003). 

Recently, observations by De Buizer (2004) and Dc 
Buizer ct al. (2003) have identified four candidate HM- 
POs that have been observed in the mid-infrared with 
high-angular resolution. These HMPO candidates lie in 
areas of high mass star formation, as evidenced by nearby 
UC H II regions. In fact in two of these four cases, high 
spatial resolution observations from 8m telescopes were 
needed to simply isolate the mid-infrared emission of the 
HMPOs from the nearby UC H II regions. This means 
that, though more sensitive than what is achievable from 
the ground, data from infrared satellites like MSX (^5" at 
8 /Ltm, ~13" at 21 /xm), and even the Spitzer Space Tele- 
scope (~2" at 8 /Ltm, ~6" at 24 /xm), are in general not 
adequate for these studies because of their relatively poor 
angular resolutions. By being able to isolate emission from 
the HMPO alone, mid-infrared observations from ground- 
based 8 to lOm-class telescopes represent a invaluable tool 
in the study of these earliest stages of massive star forma- 
tion. 

In this article we present new sub-arcsccond mid- 
infrared observations of three of the four candidate HM- 
POs identified by De Buizer (2004) and De Buizer et al. 
(2003). The fourth HMPO candidate, G305.20-F0.21, will 
be the subject of a following article. Using the Gemini 
South telescope we imaged these three HMPO candidates 
through up to 11 filters, sampling their SEDs over the 
entire 2-25 fim infrared atmospheric window. Our mo- 
tivation for performing these observations was to sample 
the SEDs of HMPOs with sufficient spectral resolution and 
coverage to accurately test new SED models. These mod- 
els use envelopes that take into account deviations from 
spherical geometry in their inner region due to rotation 



and at large scales due to natural elongation of the source. 
By matching these models to the data we attempt to "back 
out" of the models detailed physical parameters for these 
earliest stages of massive stellar life. In this way, we hope 
to be able to learn some details about the massive star 
formation process during these highly obscured stages, as 
well as the strength and limitations of the non-spherical 
SED models. 

2. OBSERVATIONS 

Observations were carried out at Gemini South over the 
time period between 2003 September and 2004 November. 
Imaging was performed with the Thermal-Region Camera 
and Spectrograph (T-ReCS). The instrument employs a 
Raytheon 320x240 pixel Si: As IBC array which is opti- 
mized for use in the 7-26 fxm wavelength range, but can 
perform with modest sensitivity in the 2-5 /xm wavelength 
range as well. The pixel scale is 0.089"/pixel, yielding a 
field of view of 28'.'8x21'.'6. Sky and telescope subtraction 
were achieved through the standard chop-nod technique. 

Images were performed with K, L, and M filters, as well 
as six narrowband silicate filters encompassing the full 10 
/im A'^-band spectral window, and two narrowband filters 
within the 20 /xm Q-band spectral window. Effective cen- 
tral wavelengths for these filters and their spectral range 
are given in Table 1. For Gll.94-0.62 and G45.07-I-0.13 
wc also have images from Dc Buizer ct al. (2003) through 
the Q3 (Ao=20.81 ^^m, AA=19.99-21.64 /xm) filter from 
the NASA Infrared Telescope Facility (IRTF) using the 
instrument MIRLIN. 

Table 1 also lists the on-source integration times through 
each of the T-RcCS filters for each HMPO candidate field. 
In addition, this table lists the standard star used to fiux 
calibrate each image, and the assumed flux density for that 
standard through the given filter. These assumed standard 
star flux densities were found by convolving the spectral 
irradiance templates of the stars from Cohen et al. (1999) 
with the given T-ReCS filter transmission profile. 

Flux densities were derived for each filter for the HMPO 
candidate as well as any other mid-infrared sources on the 
field. Tables 2-4 list the fiux densities for all sources. These 
flux densities are quoted with their l-cr total error, which 
is a quadrature addition of the statistical variation from 
the aperture photometry (due to the standard deviation of 
the background array noise) and the fiux calibration error. 
Usually the flux calibration error is found from the vari- 
ation of the standard star flux throughout the course of 
the night, however our observations were taken piece- wise 
on many nights throughout several semesters. Flux cal- 
ibration was therefore achieved by observing a standard 
star at a similar airmass just before or after the obser- 
vation of the target science field. Therefore, we do not 
have calibration variability statistics for each filter on each 
night the data were taken. However, the calibration factor 
[ratio of accepted fiux in Janskys to analog-to-digital con- 
verter units (ADUs) per second per pixel] derived from the 
standard star observations through each filter varied little 
throughout 14 months of data collection, and so we use 
here the standard deviation of the calibration factor over 
this 14 month period as the fiux calibration error. This 
can be considered an extremely conservative estimate of 
the errors, taking into account a wide variety of atmo- 
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spheric and observing conditions that affect flux calibra- 
tion. Most filters have modest flux calibration errors (for 
mid-infrared observations) with standard deviations be- 
tween 2 and 10%. Flux calibration through certain filters 
is more difficult due to the presence of various atmospheric 
absorption lines contaminating the filter bandpass, some 
of which can be highly variable. Those filters most affected 
are the 7.7 /xm (21%), 12.3 /xm (19%), 18.3 yum (15%), and 
24.6 iim (23%) filters. 

Furthermore, the continuum observed through these fil- 
ters may be affected by the presence of certain spectral fea- 
tures common to young stellar objects. In particular, the 
12.3 fim continuum flux densities for the targets may be 
skewed to higher values due to the presence of bright [Nell] 
lines at 12.8 /xm (see for instance, Faison et al. 1998). The 
18.3 fxm and 20.8 fxm observations can also be skewed high 
or low depending on the presence in either emission or ab- 
sorption of the 18 /xm silicate feature. 

3. THE TARGETS 

The three sources below were chosen because they arc 
among the only HMPO candidates known to be bright in 
the mid-infrared (De Buizer 2004, De Buizer et al. 2003) 
and are sufficiently far enough away from nearby H II re- 
gions to be studied in detail without confusion. All three 
sources are thought to be HMPOs primarily because they 
are sources found in regions of known high mass star for- 
mation yet have no UC H II emission of their own, and 
all are coincident with molecular maser emission which is 
thought to trace massive star formation and its processes. 
Each source has further unique evidence indicating that it 
is a HMPO, which we discuss in detail below. A summary 
of the observations for each object will also be discussed 
in this context. 

3.1. Gil. 94-0.62 

Of the four sources identified as mid-infrared bright 
HMPO candidates by De Buizer (2004) and De Buizer et 
al. (2003), the status of the source in G11.94-0.62:DRT03 
1 as a genuine HMPO is the most uncertain. However 
the HMPO candidate in Gil. 94-0. 62 satisfies, at the very 
least, the; minimum criteria outlined at the beginning of 
this section. DRT03 1 was found as an unresolved mid- 
infrared source in the survey of De Buizer et al. (2003) lo- 
cated at the site of a water maser clump offset ^10" from 
a cometary UC H II region. DRT03 1 has no detected cm 
radio continuum emission of its own. 

However, Gil. 94-0. 62 has not been observed with high 
spatial resolution molecular line imaging. Some molecular 
tracers have been discovered toward the region, like CS 
(Bronfman et al. 1996) and NH3 (Cesaroni, Walmsley, & 
Churchwell 1992), but it is not known if this emission is 
related to the nearby UC H II region, the HMPO candi- 
date, or both. While it is encouraging that some molecules 
have been foimd in the region, CH3CN (a tracer of denser 
and hotter molecular material and indicator of the pres- 
ence of HMCs) has been searched for in the region and was 
not detected (Watt & Mundy 1999). The UC H II region 
was also imaged at 2.7 mm by Watt & Mundy (1999), but 
there was no detection at the site of DRT03 1. 

De Buizer et al. (2003) derive a lower limit to the bolo- 
metric luminosity for DRT03 1, finding that it is at least as 



bright as a B9 star. Though this spectral type is derived 
from an extreme lower limit to the bolometric luminosity, 
it does not appear that this source is a massive (M>1OM0) 
protostar, but more likely an intermediate mass protostel- 
lar object. If this is true, given its observational similari- 
ties to its higher mass counterparts, it will be interesting 
to study the similarities and differences in detail to learn 
more about the interrelation of the intermediate and high 
mass star formation processes. It would also have interest- 
ing c;onsequences for the hypotheses regarding excitation 
of water masers by lower mass (non-ionizing) stars (Forster 
& Caswell 2000). 

Our sensitive Gemini mid-infrared observations of the 
region around Gil. 94-0. 62 have revealed extended dust 
structure, the likes of which are not seen at any other spec- 
tral regime (Figure 1). These images at the sub-arcsecond 
angular resolution of Gemini confirm the results of the 
lower resolution observations of De Buizer et al. (2003) 
that the mid-infrared morphology of the UC H II region 
is not cometary, as is seen in cm radio continuum emis- 
sion (Wood & Churchwell 1989). The HMPO candidate, 
DRT03 1, is clearly detected at shorter mid-infrared wave- 
lengths but is only a 3-a and 4-a detection at 20.8 and 
24.6 /Ltm, respectively. 

With no information on the Rayleigh-Jeans side of the 
SED for the HMPO candidate alone, we will not be able 
to constrain the models at the longer wavelengths. Table 
2 lists the observed flux densities for the HMPO candidate 
as well as for the other mid-infrared sources in the field. 
It is not clear if each of these mid-infrared sources house 
a stellar component or if they are simply knots of dust in 
an extensive H II region. 

3.2. G29.96-0.02 

The HMPO candidate in G29.96-0.02 is sufficiently 
well-studied to be considered a prototype. The hot molec- 
ular component of this source was first discovered in the 
observations of Cesaroni et al. (1994). In this work, a 
hot {Tk> 50 K) core of molecular material was detected, 
traced by NH3(4,4) line emission near to. but offset ^2" 
from a cometary-shaped H II region. This hot molecu- 
lar core was found to be coincident with a cluster of wa- 
ter masers. Cesaroni et al. (1994) conjectured that this 
source, like the other HMCs in their survey, is internally 
heated by a developing massive protostar. However, as 
was mentioned in the introduction, some HMCs could be 
externally heated by the nearby massive stars powering the 
UC H II regions. In the case of G29.96-0.02, the millimeter 
line study of Gibb, Wyrowski, & Mundy (2003) provided 
strong evidence that the HMC is internally heated. There- 
fore, the HMC in G29.96-0.02 appears to be a genuine 
HMPO. 

The HMPO in G29.96-0.02 was directly imaged for the 
first time in the mid-infrared by De Buizer et al. (2002) 
on the Gemini North 8-m telescope. The HMPO lies -3" 
away from the radio continuum peak of the bright H II re- 
gion and the emission from the UC H II region peak cannot 
be resolved from the HMPO on 4-m class telescopes (Wat- 
son et al. 2003). However, the large aperture of Gemini 
has sufficient resolving power in the mid-infrared (~0'.'5) 
to pull in the peak of the dust emission from the H II re- 
gion to allow direct observation of the HMPO. Even so, 
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the resolved thermal dust emission from the H II region 
is still extensive, and therefore we cannot avoid observing 
the HMPO overlaid on a diffuse background of emission 
(Figure 2). Following the technique outlined by De Buizer 
et al. (2002), one can fit this background emission and 
subtract it off, effectively leaving only the mid-infrared 
emission of the HMPO. This technique was employed for 
the images obtained in the work presented here. In Table 
3 we present the flux densities for the HMPO candidates 
and the UC H II, as well as the flux density of the HMPO 
candidate without subtracting off the background due to 
the extended UC H II emission. This was done so that it 
is transparent to the reader how much flux has been sub- 
tracted off, and what the absolute upper limits are for the 
flux densities of the HMPO candidate. 

The G29.96-0.02 HMPO has been studied in great de- 
tail in molecular line studies (e.g., Cesaroni et al. 1998, 
Wyrowski, Schilkc, & Walmslcy 1999; Maxia et al. 2001.) 
However, other than the mid-infrared observations of De 
Buizer et al. (2002) and those presented here, only a few 
continuum observations at other wavelengths exist. This 
source has been observed at 1 and 3 mm by Maxia et 
al. (2001). However, as pointed out by De Buizer et al. 
(2002), there appears to be problems with these flux den- 
sities since they lead to a negative value (see Figure 6) for 
the index of dust opacity, /?, whereas most models assume 
l</3<2 for A>200 /zm. 

Another observation at 3 mm was carried out by Olmi cl 
al. (2003). They did not resolve the HMPO, and therefore 
there might be possible contamination from the nearby 
UC H II region. They try to obtain the flux density of 
the HMPO at 3 mm, making a subtraction of the free-free 
contamination estimated from the 2 cm image. This mea- 
surement together with the detection of the HMPO at 1.4 
mm carried out by Wyrowski et al. (2005, personal com- 
munication) indicate that G29.96-0.02 HMPO is a strong 
emitter at millimeter wavelengths. 

3.3. G45. 07+0.13 

Listed as a "known hot core" in the review article of 
Kurtz et al. (2000), the HMPO in G45.07-F0.13 gained 
attention mostly through the work of Hunter, Phillips, 
& Menten (1997). However most of the observations by 
Hunter, Phillips, & Menten (1997) are low spatial reso- 
lution continuum and molecular line maps observed to- 
wards UC H II regions. In the case of G45.07-I-0.13 they 
found several molecular tracers, some believed to be in out- 
flow, emanating from the location of the UC H II region. 
Hunter, Phillips, & Menten (1997) contest that there is a 
single star here, embedded in a molecular core, at a very 
early stage of development such that it has just begun to 
develop a UC H II region. For these reasons it has been 
labeled a hot core. 

The observations of De Buizer et al. (2003) showed that 
there are three mid-infrared sources at this location, all 
within 8" of each other. The brightest mid-infrared source 
on the field is coincident with the UC H II region. How- 
ever, at a location ~2" north of the UC H II region there 
is a group of water masers associated with a mid-infrared 
source that has no cm radio continuum emission (Figure 
3). De Buizer et al. (2003) showed that this northern mid- 

^ The coordinates for the HMPO candidate are wrong in De Buizer et 



infrared source, DRT03 3, has a very bright mid-infrared 
luminosity of ~4x 10^ Lq. Assuming that the mid- infrared 
luminosity is equal to the bolometric luminosity (a gross 
underestimate), DRT03 3 has the equivalent luminosity of 
a BO star. We therefore have a source luminous enough to 
be a massive protostar, situated at the location of a group 
of water masers, yet this source has no cm radio contin- 
uum emission of its own. For these reasons, De Buizer et 
al. (2003) claim that DRT03 3 is the true HMPO in this 
ficldi. 

Although there are multiple sources of thermal dust 
emission present here, it does appear from the high spa- 
tial resolution CS maps of Hunter, Phillips, and Menten 
(1997) that the stellar source powering the UC H II re- 
gion is the only viable source of the outflow in the region 
(contrary to what was stated by De Buizer et al. 2003). 
Since no such tracers are seen coming from the HMPO 
candidate, this perhaps confirms the youth of the object, 
placing it in a stage of formation before the onset of an 
outflow. While the resolution of the 800 and 450 /zm maps 
of Hunter, Phillips, & Menten (1997) are too coarse to as- 
certain which mid-infrared source may be dominating the 
emission at these longer wavelengths, their higher resolu- 
tion (~2") 3 mm map shows no emission from the location 
of DRT03 3, though emission is detected towards the other 
two mid-infrared sources on the fleld. 

Our new observations of this field with Gemini have re- 
vealed the presence of four more faint mid-infrared sources 
present (Figure 3) in addition to those already found by De 
Buizer et al. (2003). It appears that this may be a small 
cluster of star formation centered on the 09 star powering 
the UC H II region. DRT03 3 is seen at all 11 mid-infrared 
wavelengths, as is the UC H II region (DRT03 2). Table 4 
lists the flux densities for all mid-infrared sources on the 
field. Unfortunately, we have no detections of the HMPO 
candidate alone at any wavelength on the Rayleigh- Jeans 
side of the SED, and have to rely upon integrated flux den- 
sity measurements from the whole region (Su et al. 2004) 
as upper limits to constrain the models. 

4. MODELING 

4.1. Description of the Model 

We will calculate the SED of a HMPO by modeling it 
as an envelope of gas and dust that is freely falling onto a 
recently formed massive central star, which is responsible 
for the heating. Osorio, Lizano & D'Alessio (1999) mod- 
eled the SED of a sample of such objects by assuming that 
the envelope has spherical symmetry. For that work, only 
upper limits were available for the mid-infrared flux densi- 
ties of most of the sources. Because this wavelength range 
is very sensitive to the geometry of the source, here we ex- 
tend that work by including rotation and flattening of the 
envelope in the models, so that the results can be tested 
with the new mid-infrared data presented in this paper. 
For this purpose, we have followed procedures similar to 
those already developed in low-mass star models (Osorio 
et al. 2003). 

Envelopes with a density distribution such as that given 

by Terebey, Shu, & Cassen (1984, hereafter TSC en- 
velopes) assume that the rotation of the infalling material 

al. (2003) . See the caption of Figure 3 for correct coordinates 
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becomes important only in the vicinity of the centrifugal 
radius, Rc = TqI^o/GM*, where f^o is the angular velocity 
(assumed to be constant and small) at a distant reference 
radius ro, and M* is the mass of the central protostar. 

The TSC envelopes arc flattened only in the inner re- 
gion, being essentially spherical in the outer region. The 
inner regions (r < ro) of these envelopes are described 
by the solution of Cassen & Moosman (1981) and Ulrich 
(1976) (hereafter referred as CMU). The solution is de- 
termined by specifying the velocities and densities at the 
(distant) reference radius ro, where the velocity is nearly 
equal to the radial free-fall value, with only a small az- 
imuthal component corresponding to a constant angular 
velocity, and the density is set by the required mass infall 
rate M. With these assumptions, the density of the inner 
envelope (r < ro) is given by 

^°^"^''^) = 4.(GAf,r3)i/2 

\ cos 6o J \COS^o / ' 

where 6 and Oq are the position angles with respect to 
the rotational symmetry axis of the system at radial dis- 
tances from the central mass r and ro, respectively. At 
large distances (r » Rc) the motions are radial {9 Oq) 
and the density tends to the density distribution for ra- 
dial free-fall at constant mass infall rate. For r < i?c, the 
motions become significantly non-radial, and the angular 
momentum of the infalling material causes it to land on 
a disk {9 = 90°) at distances < r < Rc. Therefore, Rc 
is the largest radius on the equatorial plane that receives 
the infalling material. 

A more complex description can be given in terms of 
intrinsically flattened envelopes, such as those with the 
density distribution resulting from the gravitational col- 
lapse of a sheet initially in hydrostatic equilibrium (Hart- 
mann et al. 1994; Hartmann, Calvet, & Boss 1996). These 
envelopes are flattened not only in their inner region, be- 
cause of rotation (as the TSC envelopes), but also at large 
scales due to the natural elongation of the cloud. If the 
rotational axis is perpendicular to the plane of the self- 
gravitating layer, the density distribution of the infalling 
material will be axisymmetric and can be written as: 

p{r,9) = pcMv{r,9)r]cosh.~'^{r]cos9o)tanh~^{r]), 

where pcMU is the CMU density distribution for the same 
M and Rc, and 77 is a measure of the degree of flatten- 
ing (Hartmann et al. 1996). The degree of flattening in 
these envelopes is measured by the parameter rj = Ront/H, 
where i?out is the outer radius of the envelope and H is 
the scale height. Hereafter, we will designate them as r]- 
envelopes. 

The scale of this density distribution can be set by using 
the parameter pi, defined as 

M 

" 47r(2GM,r3)i/2' 

where r\ is a reference radius. This reference density cor- 
responds to the density that a spherically symmetric free- 
falling envelope with the same mass accretion rate, M, 
would have at the reference radius ri. In our modeling, 



we will use ri = 1 AU, and we will designate the reference 
density as pi au • 

The remaining parameters that define the density dis- 
tribution of the ?7-envelopes are the centrifugal radius, Rc, 
and the inclination of the polar (rotational) axis to the 
line of sight, i. Since ??-envelopes are fiattened not only 
in their inner region but also at large scales most of their 
material is accumulated over the equatorial plane, having 
an extinction lower than a TSC envelope of the same mass 
provided the line of sight is not oriented too close to the 
edge-on position. Therefore, a flattened ?7-envelope model 
predicts stronger near and mid-infrared emission and less 
deep absorption features than those predicted by a TSC 
envelope model of the same mass. Given that our sources 
are strong mid-infrared emitters, we will preferentially use 
the ry-envelope density distribution to model our sources. 

The temperature distribution along the envelope, as well 
as the SED from 2 /im to 2 cm, arc determined following 
the procedures described by Kenyon, Calvet, & Hartmann 
(1993, hereafter KCH93), Calvet et al. (1994), Hartmann 
et al. (1996), and Osorio et al. (2003). For a given lumi- 
nosity of the central star, L*, which is assumed to be the 
only heating source (i.e., neglecting other possible sources 
of heating such as accretion energy from a disk or the enve- 
lope), the temperature structure in the envelope is calcu- 
lated from the condition of radiative equilibrium, using the 
approximation of an angle-averaged density distribution. 
Nevertheless, the c!mc^rgc!nt flux density is calculated rising 
the exact density distribution. Because we want to com- 
pare our results with high spectral resolution mid-infrared 
data, our model spectra wavelengths were finely sampled 
(AlogA — 0.02) in the range 3-25 /xm in order to define 
with high spectral resolution the 10 fiuv silicate absorption 
feature. Outside this range, wavelengths were sampled at 
AlogA = 0.2-0.3. Additionally, in order to make a proper 
comparison with the observations, we convolved our re- 
sults with the T-ReCS filter transmission profile. Given 
that the spectral response of the filters is fairly constant, 
we do not find large variations (only ~5%) between the 
averaged and original spectrum. 

In our models we use a recently improved dust opacity 
law, corresponding to a mixture of compounds that was 
adjusted by matching the observations of the well known 
low-mass protostar L1551 IRS5 (Osorio et al. 2003). This 
mixture includes graphite, astronomical silicates, troilite 
and water ice, with the standard grain-size distribution of 
the interstellar medium, n{a) oc a~^^^, with a minimum 
size of 0.005 fim and a maximum size of 0.3 /im. The 
details of the dust opacity calculations are described by 
D'Alessio, Calvet, & Hartmann (2001). 

The inner radius of the envelope, is taken to be 
the dust destruction radius, that is assumed to occur at a 
temperature of 1200 K, corresponding to the sublimation 
temperature of silicates at low densities (D'Alessio 1996). 
The centrifugal radius is assumed to be considerably larger 
than the dust destruction radius in order to have a signifi- 
cant degree of flattening due to rotation in the inner enve- 
lope. The outer radius of the envelope, i?out, is obtained 
from our 11.7 /im images. 

In summary, once the dust opacity law and the inner 
and outer radius of the envelope are defined, the free pa- 
rameters of our model are r}, piAu, Rc, i, and L*. In 
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the following, we will discuss on the contribution of these 
parameters to the resulting SED. 

4.2. Behavior of the SED 

In order to illustrate the effect of the flattening of the 
envelope, we have calculated the resulting SED for a TSC 
envelope and compared it with an jj = 2 envelope of the 
same mass. The results are shown in Figure 4, where it can 
be seen that the flattened 77 = 2 envelope predicts stronger 
near and mid-infrared emission, and less deep absorption 
features than the TSC envelope, both for low {i = 30°) and 
moderately high {i = 60°) inclination angles. This behav- 
ior is a consequence of the flatter distribution of material 
in the r] = 2 envelope, resulting in an extinction lower 
than the TSC envelope over a wide range of viewing an- 
gles (except for a nearly edge-on view, i = 90°, where the 
behavior would be the opposite). Thus, in general, the 
flattened ?7-envelopes allow the escape of a larger amount 
of near and mid-infrared radiation, and consequently they 
predict stronger flux densities than the TSC envelopes in 
this wavelength range. Since our sources are strong in- 
frared emitters, in our models we will use ?7-cnvelopes with 
a significant degree of flattening, 77, of the order of 2. 

The observed SEDs of our sources are well sampled in 
the near and mid-infrared wavelength range. However, this 
wavelength coverage in insufficient for the models to deter- 
mine the value of the bolomctric luminosity (that we take 
to coincide with the stellar luminosity, L,). The bolomct- 
ric luminosity is well constrained by the far-infrared data 
(where the peak of the SED occurs), and to a less extent 
by the submillimeter or millimeter data. Unfortunately, 
the available far-infrared data for our sources do not have 
the appropriate angular resolution to avoid contamination 
from other nearby sources, and should be considered as 
upper limits, while millimeter data are only available for 
one of our sources. Therefore, in general, we constrain the 
value of the luminosity by fitting the mid-infrared data 
closer to the SED peak, that are helpful to set a lower 
limit for the luminosity. 

On the other hand, the depth and shape of the sili- 
cate absorption feature (which is very well delineated by 
our mid-infrared data) determines the extinction along the 
line-of-sight, that in turn is related to the inclination an- 
gle, i, the density scale, pi, and the centrifugal radius. Re- 
The depth of the 10 /im silicate absorption increases not 
only with the density scale of the envelope, but also with 
the inclination angle, since the extinction is much smaller 
along the rotational pole (i = 0) than along the equator 
{i = 90°) (see Osorio et al. 2003). Thus, the behavior 
of the parameters pi and i is somewhat interchangeable 
in the mid-infrared (although such behavior is different at 
millimeter wavelengths) . For instance, an envelope with a 
combination of low density and high inclination could have 
an extinction and mid-infrared emission similar to an enve- 
lope with a higher density but viewed pole-on. Therefore, 
for the sources where millimeter data are not available, it 
is difficult to constrain i and pi simultaneously, and we will 
search for fits both for high as well as for low inclination 
angles. 

The behavior of the SED with the centrifugal radius, i?c, 
is more complex. Increasing the value of Rc will widen the 
region of evacuated material, producing an overall decrease 



of the extinction along the line of sight, but, at the same 
time, it will result in a decrease of the temperature of the 
envelope, because the material tends to pile up at larger 
distances from the source of luminosity. Thus, depending 
on the wavelength range, one of the two effects will be the 
dominant. At short wavelengths, where the dust opacity 
is higher, it is expected that the lower extinction of the 
envelopes with large values of Rc will result in an increase 
of the observed emission. At longer wavelengths, where 
the opacity is low enough (except in the silicate absorp- 
tion feature), the dominant effect is expected to be due 
to the decrease in the temperature, resulting in a decrease 
of the observed emission for larger values of Rc- This is 
illustrated in Figure 4, where it is shown that the emission 
of the envelope with Rc = 600 AU dominates in the near- 
infrared range, while the emission of the Rc = 250 AU en- 
velope dominates at longer wavelengths. As shown in the 
figure, this result is valid both for low and relatively high 
inclination angles, although we note that the point where 
the small Rc envelope becomes dominant shifts slightly 
to longer wavelengths when increasing the inclination an- 
gle (from ~ 2 /xm for i = 30° to ~ 5 /zm for i = 60°). 
This shift occurs because envelopes seen at higher incli- 
nation angles are more opaque, thus requiring a longer 
wavelength (where the opacity is lower) for the small Rc 
envelope to become dominant. 

Since we expect the value of Rc to be of the order of 
the radius of a possible circmnstcllar disk, in our mod- 
eling we will explore a range of values for Rc from tens 
to hundreds of AUs, depending on the luminosity of the 
source. The upper limit of this range is suggested by the 
highest angular resolution observations of disks reported 
towards massive protostars (e.g., Shepherd, Claussen, & 
Kurtz 2001), while the lower limit is set by the smallest 
values recently reported for protoplanetary disks around 
low- mass YSOs (Rodriguez et al. 2005). 

From the values of the parameters obtained from our 
fits to the observed SEDs we can derive other physical pa- 
rameters, such as the total mass of the envelope, M^nv, 
obtained by integration of the density distribution, and 
the mass of the central star, M^, derived from its lumi- 
nosity, L*. Once the mass of the star is known, the mass 
accretion rate in the envelope, M, can be derived from the 
value of pi. 

Yorke (1984) and Walmsley (1995) proposed that some 
objects could be accreting material at a rate high enough 
to quench the development of an incipient HII region. 
These objects would lack detectable free-free emission, be- 
ing good candidates to be the precursors of the UCHII re- 
gions. In the case of spherical accretion, it is easy to show 
that the critical mass accretion rate required to reduce the 
HII region to a small volume close to the stellar surface is 
-^-^crit = (STrGmfj^a^^ NiM^y^^ , where mn is the hydrogen 
mass, and a is the recombination coefficient (excluding 
captures to the n = 1 level). We will adopt the value, 
Merit, as a critical measure of the mass accretion rate re- 
quired to "choke off" the possible development of an UC 
H II region. Since we are assuming departures from spher- 
ical symmetry in our envelopes, in reality a more complex 
calculation would be required to accurately derive the cor- 
responding critical mass rate. However, for simplicity we 
will assume that if the ratio between the mass accretion 
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rate in our envelope and the critical spherical mass ac- 
cretion rate M/Mcrit 1, it would indicate that pho- 
toionized emission is not expected to be detectable in the 
source. 

5. COMPARISON OF THE MODELS WITH THE 
OBSERVATIONS 

In this section, we estimate the physical parameters of 
the three HMPO candidates by fitting our models to the 
observed SEDs, using also the available information on 
size and morphology inferred from the images. Given that 
spherical envelopes do not appear to be able to account for 
the observed properties of the SEDs in the mid-infrared, 
wc have modeled the soiirccs using r^-cnvclopcs with a sig- 
nificant degree of elongation (see discussion in §4.2). Since 
the temperature distribution in the envelopes is estimated 
using the approximation of an angle- averaged density dis- 
tribution, valid for moderately elongated envelopes, we did 
not attempt to consider very high values of ry, that corre- 
spond to extremely elongated envelopes where the temper- 
ature calculations probably would not be accurate enough. 
Therefore, we explored values of rj around 2, and in all 
three sources wc were able to get good fits for rj — 2.5. 
We note that this value coincides with the one that pro- 
vided the best fit for the low-mass protostar L1551-IRS5, 
after exploring a wide range of values of the rj parameter 
(Osorio ct al. 2003). 

Given the incomplete and non-uniform coverage of the 
observed SED, the fitting process cannot be automated 
and should be carried out manually on a case-by-case ba- 
sis. As was already mentioned in the previous section, the 
luminosity is well constrained by the peak of the SED, that 
is expected to happen in the far-infrared. Unfortunately, 
at present, far- infrared observations do not reach the angu- 
lar resolution required to properly isolate the emission of 
these distant sources from contamination of other sources 
in their vicinity. Therefore, in general, the far-infrared 
data are considered only as upper limits and we have to 
use the mid-infrared data to obtain a range of possible 
values for the luminosity. Thus, the general strategy was 
to first run a set of models with increasing luminosities 
until we found a value that was able to reproduce the ob- 
served mid-infrared flux density near the peak of the SED 
and consistent with the millimeter and submillimeter data. 
Once an approximate value of the luminosity was found, 
we proceeded to run a set of models to find the density 
scale, which is then basically determined by the millime- 
ter and submillimeter data points. When the density scale 
was fixed, we ran models with different inclinations in or- 
der to find the best fit to the absorption feature in the 
mid-infrared, which is particularly sensitive to the value 
of the inclination angle. Finally, we refined the model- 
ing by testing a range of possible values of the centrifugal 
radius. 

If the data point coverage of the SED is good enough 
we should end this process with essentially one final model 
(the "best fit model") that is able to reproduce the ob- 
served data points within the constrains imposed by the 
observational uncertainties. In the case that only upper 
limits are available in the millimeter and submillimeter 
domain, the density scale and inclination angle cannot be 
constrained simultaneously, and in this case we present 



two fits, one for a low inclination angle, and a second one 
for a high inclination angle. The value of the centrifugal 
radius not well constrained in this case either. However, 
on the basis of an additional analysis of the quality of the 
data points and the peculiarities of each source (see dis- 
cussion of individual sources) we attempt to favor one of 
the fits as the best fit model. 

Table 5 lists the values of the parameters of our best fit 
models: the inner, outer and centrifugal radii of the enve- 
lope, the stellar luminosity, the refc^rence density at 1 AU, 
and the inclination angle. Our best-fit model SEDs are 
shown in Figures 5, 6, and 7. Table 6 lists other derived 
parameters for each source: the mass of the envelope, its 
temperature at a radius of 1000 AU, the mass and spectral 
type of the central protostar, the mass accretion rate, the 
rate of ionizing photons, and the critical mass accretion 
rate. These parameters have been derived following the 
procedures explained in the footnotes to the table. 

5.1. Gil. 94-0.62 HMPO 

This HMPO has been observed through continuum mea- 
surements at millimeter (Watt & Mundy 1999), submil- 
limeter (Walsh et al. 2003), and mid- infrared wavelengths. 
The submillimeter data will be used in our modeling only 
as upper limits, since they have been obtained with a large 
beam (10"-18") and therefore it is not possible to separate 
the emission of the nearby UC H II region from that of the 
HMPO. Likewise, the millimeter fiux densities will be used 
also as upper limits, since the analysis of Watt & Mmidy 
(1999) suggests that the millimeter emission observed to- 
wards this source corresponds to free-free emission of the 
nearby UC H II region with a negligible contribution from 
dust. 

In general, to constrain the physical parameters of the 
sources, we use the observed SED as well as the informa- 
tion on size and morphology given by the images. How- 
ever, the Gil. 94— 0.62 HMPO source appears unresolved 
in the near and mid- infrared images (see Figure 1) and 
therefore, we have only an upper limit for its size. 

Because we considered the far-infrared, siibmillimeter 
and millimeter data only as upper limits, the luminos- 
ity of this source is constrained essentially by the data 
at wavelengths near 20 /xm, that indicate that this is a 
low luminosity source. We obtain a value of 75 Lq for the 
luminosity of this object, corresponding to a AO star of 5 
Mq. Since we only have upper limits for the flux den- 
sities in the submillimctcr-millimeter wavelength range, 
we cannot constrain the density and inclination simul- 
taneously. Therefore, we searched for fits both at low 
(?;=20°-30°) and high inchnation angles (i=50°-60°). At 
low inclination angles, we obtained a quite good fit for 
i =30°, piAu = 7.5 X 10"" g cm-3, and i?c = 30 AU 
(dot-dashed line in Figure 5), but this model predicts too 
much emission at 18.3 and 24.6 /xm. In order to repro- 
duce the observational data in the 18-25 /im wavelength 
range, we require models with a high inclination angle. 
Our best fit model is obtained for a high inclination angle, 
i = 53°, with piAu = 2x10-13 g cm-3 (Me„„ = 1 M©), 
and Rc = 30 AU (solid line in Figure 5). We can ob- 
tain reasonable fits for 1 x 10^^^ < piAU < 4 x 10-^^ g 
cm-3 and 20 < i?c < 200 AU. For example, a model with 
Rc = 100 AU is also feasible (dotted line in Figure 5). 
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However, we favor the Rc = 30 AU model (solid line in 
Figure 5, and Tables 5 and 6) because it reproduces bet- 
ter the flux densities shortward of 8.7 /iin and the depth 
of the 10 /xm silicate absorption feature. 

For our favored model, the ratio M / Mcnt ^ 1 (sec 
Table 6). Therefore, the source could not develop a de- 
tectable H II region. Furthermore, since the luminos- 
ity obtained in this model is quite low, free-free emis- 
sion from the HMPO would be hardly detectable, even 
in the case of a negligible accretion rate. Thus, Gil. 94- 
0.62 is not luminous enough to be considered as a HMPO. 
Most likely, given the low temperature of the envelope 
(TioooAU = 40 K), it would be better considered as a 
"warm core" that is forming an intermediate-mass star, 
with an infall rate more than one order of magnitude 
higher than the typical infall rates of low- mass protostars 
(see KCH93). 

5.2. G29.96-0.02 HMPO 

The G29.96-0.02 HMPO is the only source of our sam- 
ple that appears to be associated with strong millimeter 
emission. Thus, for this source we have both infrared and 
millimeter data to constrain better the models. The strong 
millimeter omission siiggests a dense or luminous envelope. 
In addition, the relatively strong near-infrared emission 
observed suggests that the envelope should be seen nearly 
pole-on, since in this direction there is a decrease of the 
density with respect to the equatorial plane that allows 
the escape of near-infrared radiation. 

A reasonable fit (solid-line in Figure 6), explaining si- 
multaneously the strong millimeter emission as well as 
the near-infrared flux densities, is obtained by assuming 
an extremely dense (piAU = 3.1 x 10~^^ g cm~^), and 
therefore very massive envelope {Menv = 576 Mq), with 
a centrifugal radius of Rc = 570 AU, seen with a low 
inclination angle {i = 12°), and heated by a luminous 
(L* = 1.8 X 10'' Lq) B1 star (see Tables 5 and 6). This 
value of the luminosity is a lower limit, since it is the small- 
est value that is consistent with the mid-infrared and mil- 
limeter data points, and we estimate that it can be higher 
up to a factor of three. Likewise, we estimate that the un- 
certainty in the density scale is ~^30%, and that the range 
of plausible values is 10° < i < 20° for the inclination an- 
gle, and 500 < Rc < 1200 AU for the centrifugal radius. 
Note that the 9.7 /im data point is only an upper limit, 
and therefore it does not set the depth of the silicate ab- 
sorption feature. Also, we want to point out that we have 
not attempted that our models reproduce the exact values 
of the observational data points in the millimeter range, 
due to the problems discussed in §3.2 and because errors 
are not available for some of these data points. 

Our best fit model predicts that the flux density at 
3.9 jum should be lower than the value observed. We 
think that the excess of flux density observed at wave- 
lengths shorter than 3.9 /im could be due to scattered 
light in a cavity carved out by the outflow associated with 
G29.96-0.02 HMPO (Gibb et al. 2004). In fact, Gibb 
et al. (2004) suggest that scattered fight in G29.96-0.02 
HMPO could be dominant even at 18 /xm. However, our 
models show that the emission from 4-18 fim can be re- 
produced as dust emission from an infalling flattened en- 
velope, and that scattered light may only be dominant for 



wavelengths shorter than ~4 /xm. 

From Table 6 we see that the ratio M/ Merit = 22, 000. 
Thus, we do not expect the development of an H II re- 
gion in the G29 HMPO. This is consistent with the lack 
of strong free- free emission at the position of the source. 

In order to illustrate the relevance of the millimeter con- 
straints, in Figure 6 we show an additional model (dotted- 
linc), obtained by fitting only the mid- infrared observa- 
tions. For this model, we obtain a value of the luminosity 
of only 2400 Lq, a density scale of 5.3 x 10~^^ g cm""^ and 
an envelope mass of 98 Mq. This set of parameters fits 
very well the mid-infrared data; however the model pre- 
dicts far too little emission at long wavelengths, underes- 
timating by almost one order of magnitude the luminosity 
and mass. These results illustrate that both millimeter as 
well as mid-infrared data are required in order to deter- 
mine reliable values of the derived parameters. 

It is worthwhile to mention that our models can fit the 
observed SED with a value of Rc — 570 AU, considerably 
smaller than the size of the elongated structures observed 
towards the G29.96-0.02 HMPO (> 10,000 AU, Olmi et al. 
2003) and some other high-mass protostars (1000-30000 
AU, see Cesaroni 2005 and references therein). These 
structures have been interpreted by some authors as trac- 
ing large-scale disks or rotating toroids. Our results sug- 
gest that these structures could be naturally explained as 
infalling flattened envelopes (with sizes of thousands of 
AUs), rather than as large circumstellar disks. The for- 
mation of such disks is expected to occur at a scale of the 
order of Rc, which according to our models could be of the 
order of hundreds of AUs. 

5.3. G45.07+0.13 HMPO 

In order to compare the predicted SED of the 
G45. 07-1-0. 13 HMPO with the data, we use the existing 
millimeter observations (taken from Su et al. 2004) only 
as upper limits, since they do not have an angular reso- 
lution high enough to separate the HMPO emission from 
that of the UC H II region. Therefore, we will restrict our 
comparison to our mid-infrared observations. 

The G45.07-f0.13 HMPO exhibits an anomalous spec- 
trum because in spite of its deep silicate absorption (in- 
dicative of a large extinction), it presents a very strong 
emission in the near- infrared range (4-8 /xm). For this 
soiirce, we also have considered models for both high and 
low inclination angles. A model with a highly inclined en- 
velope {i = 58°), a density scale of piAu = 7.5 x 10~^^ 
g cm~^, and a centrifugal radius of Rc =270 AU (dotted 
line in Figure 7) predicts too little emission at 4.7 /xm and 
shorter wavelengths, and also presents a silicate absorp- 
tion feature that is too weak. Alternatively, an envelope 
with a lower inclination (i = 35°), a larger density scale 
(piAu = 5.3 x 10~^^ g cm~^), and a centrifugal radius of 
Rc =370 AU (solid line in Figure 7 and Table 5) reproduces 
better the emission at 4.7 /xm as well as the depth of the sil- 
icate absorption, although it still has problems in explain- 
ing the emission at the shortest wavelengths (~ 2 /xm). 
Both models require the same luminosity (L* = 2.5 x 10^ 
Lq), implying a Bl star of 13 Mq. 

We favor the low inclination model (solid line in Fig- 
ure 7 and Tables 5 and 6) because it can explain reason- 
ably well most of values of the observed flux density. As 
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we discussed for the G29. 96-0.02 HMPO case, the excess of 
flux density observed at wavelengths shorter than 3.9 /im 
could be due to scattered light. By running other models 
with slightly different values of the parameters, and tak- 
ing into account the upper limits set by the submillimeter 
and millimeter data, we estimate uncertainties of ±5° in 
the inclination angle, a factor of two in the density scale 
and the centrifugal radius, and a factor of three in the 
luminosity with respect to the values given in Table 5. 

For the values given in Table 6, a ratio M/ Merit = 3600 
is derived. Thus, as in the previous sources, the mass ac- 
cretion rate for the G45. 07-1-0. 13 HMPO is high enough to 
quench the development of an UC H II region, that is con- 
sistent with the lack of detected free-free emission towards 
this HMPO. 

6. SUMMARY 

We have obtained mid-infrared data for four HMPO 

candidates at Gemini Observatory. These data have the 
subarcsecond angular resolution necessary to isolate the 
emission of the HMPO candidates from that of nearby 
sources. The data are well-sampled across the entire 2 to 
25 /Ltm atmospheric windows, and particularly throughout 
the 10 /im absorption feature. 

Since this wavelength range is very sensitive to the 
source geometry, we have constructed a grid of SED mod- 
els for the physical conditions (density, temperature, lumi- 
nosity) expected for the early stages of massive star forma- 
tion, and taking into account for the first time the rotation 
as well as the natural flattening of the infalling envelopes. 
These models reach a degree of complexity comparable to 
those developed for low-mass protostars. Using our mid- 
infrared data, as well as the far-infrared and millimeter 
constraints available in the literature, we have fit the SEDs 
of the HMPO candidates in Gil. 94-0.62, G29.96-0.02, 
and G45.07-I-0.13, inferring physical parameters for the 
infalling envelopes as well as for the central stars. 

Our main conclusions are the following: 

1. The HMPO candidate in Gil. 94-0.62 does not 
seem to be luminous enoiigh to be considered as a 
HMPO. Most likely, it would be better considered 
as a "warm" core that is forming an intermediate- 
mass star, with an infall rate more than one order 
of magnitude higher than the typical infall rates of 
low-mass protostars. 

2. The candidates G29.96-0.02 and G45.07-h0.13 ap- 
pear to be genuine HMPOs with high luminosi- 
ties, ~ 20000 Lq, and high accretion rates, 
M ~ lO^'^-lO^^ Mq yr~^. These values of the mass 
accretion rate cixcciod by several orders of magni- 
tude the critical value required to quench the de- 
velopment of an UC H II region. Therefore, these 



sources appear to be in a very early stage, previous 
to the development of an UC H II region. 

3. Our models are able to fit the observed SEDs 

of the HMPOs with values of the centrifugal ra- 
dius of a few hundred AUs, which arc considerably 
smaller than the size of the elongated structures 
observed towards some high-mass protostars (1000- 
30000 AU). These structures have been interpreted 
by some authors as tracing large-scale disks or ro- 
tating toroids. Our results suggest that these struc- 
tures could be naturally explained as infalling fiat- 
tened envelopes (with sizes of thousands of AUs), 
rather than as large circumstellar disks. The for- 
mation of such disks is expected to occur at a scale 
of the order of i?c, which according to our models 
could be of the order of hundreds of AUs. 

4. Our detailed modeling suggests that massive star 

formation can proceed in a way very similar to the 
formation of low-mass stars. Unfortunately, mas- 
sive protostars are more distant than low-mass pro- 
tostars, and the instrumentation currently avail- 
able in the far-infrared and millimeter ranges does 
not reach the angular resolution necessary to sepa- 
rate the HMPO emission from that of other nearby 
sources. Since these wavelength ranges are very im- 
portant to constrain properly the physical parame- 
ters of the sources, such as the luminosity, in general 
our fit models are not unique. New data in these 
wavelength domains would be most valuable to de- 
termine the physical conditions in the early stages 
of massive star formation. 



Based on observations obtained at the Gemini Obser- 
vatory, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under a cooper- 
ative agreement with the NSF on behalf of the Gem- 
ini partnership: the National Science Foundation (United 
States), the Particle Physics and Astronomy Research 
Council (United Kingdom), the National Research Coun- 
cil (Canada), CONICYT (Chile), the Austrahan Research 
Council (Australia), CNPq (Brazil) and CONICET (Ar- 
gentina). Gemini program identification numbers associ- 
ated with the data presented here are GS-2003B-DD-7 and 
GS-2004A-Q-7. M. O. acknowledges support from lAU 
Peter Gruber Foundation, Junta de Andalucfa, MCYT 
(AYA2002-00376, including FEDER funds), and AECI. 
We thank G. Anglada for useful discussions and helpful 
suggestions. We also thank F. Wyrowski for communicat- 
ing us his 1.4 mm data on G29.96— 0.02 before publication. 
Finally, we would also like to thank our referee for his/her 
thoughtful review of our manuscript and useful comments. 

Facilities: Gemini South (T-ReCS). 



REFERENCES 



Behrend, R. & Macdcr. A. 2001, A&A, 373, 190 
Bronfman, L., Nymaii, L.-A., & May, J. 1996, A&AS, 115, 81 
Calvet, N., Haitmann, L., Kenyon, S. J., & Whitney, B. A. 1994, 
ApJ, 434, 330 

Cesaroni, R., Walmsley, C. M., & Churchwell, E. 1992, A&A, 256, 
618 



Cesaroni, R., Churchwell, E., Hofner, P., Walmsley, C. M., & Kurtz, 

S. 1994, A&A, 288, 903 
Cesaroni, R., Hofner, P., Walmsley, C. M., & Churchwell, E. 1998, 

A&A, 331, 709 



10 



De Buizer, Osorio, and Calvet 



Cesaroni, R. in Dense Molecular Gas around Protostars and in 

Gala<;tic Nuclei, ed. Y.Hagiwara, W. A. Baan, & H. J. van 

Langevelde, 2005, Ap&SS, 295,5 
Churchwell, E. 1991, in The Physics of Star Formation & Early 

Stellar Evolution, ed. C. J. Lada & N. D. Kylafis (Dordrecht: 

Kluwer), 221 

Cohen, M., Walker, R. G., Carter, B., Hammersley, P., Kidger, M., 

& Noguchi, K. 1999, AJ, 117, 1864 
D'Alcssio, P., Calvct, N., & Hartmann, L. 2001, ApJ, 553, 321 
D'Alessio, P. 1996, Ph.D. thesis, Univ. Nacional Autonoma de 

Mexico 

De Buizer, J. M., Watson, A. M., Radomski, J. T., Piiia, R. K., & 

Telesco, C. M. 2002, ApJ, 564, LlOl 
De Buizer, J. M., Radomski, J. T., Telesco, C. M., & Piiia, R. K. 

2003a, ApJ, 598, 1127 
De Buizer, J. M. 2004, in lAU Symp. 221, Star Formation at High 

Angular Resolution, ed. R. Jayawardhana, M. G. Burton, & T. L. 

Bourke (San Francisco; ASP), 181 
De Buizer, J. M., Radomski, J. T., Telesco, C. M., & Piiia, R. K. 

2005, ApJS, 156, 179 
De Free, C. G., Rodriguez, L. F., & Goss, W. M. 1995, Revista 

Mexicana de Astronomia y Astrofisica, 31, 39 
Doyon, R. 1990, Ph.D. thesis. Imperial College Univ. London 
Faison, M., Churchwell, E., Hofner, P., Hackwell, J., Lynch, D. K., 

& Russell, R. W. 1998, ApJ, 500, 280 
Forster, J. R., & Caswell, J. L. 2000, ApJ, 530, 371 
Gibb, A. G., Wyrowski, F., & Mundy, L. G. 2003, in SFChom 2002: 

Chemistry as a Diagnostic of Star Formation, ed. C. L. Curry & 

M. Fich (Ottawa, NRC Press), 214 
Gibb, A. G., Wyrowski, F., & Mundy L. G. 2004, ApJ, 616, 301 
Hartmann, L., Boss, A., Calvct. N.. & Whitney B. 1994, ApJ, 430, 

L49 

Hartmann, L., Calvct, N., & Boss, A. 1996, ApJ, 464, 387 
Hofner, P., & Churchwell. E. 1996, A&AS, 120, 283 
Hofner, P., Wicscmeycr, H., & Hcnning, T. 2001, ApJ, 549, 425 
Hunter, T. R., Phillips, T. G., & Mcntcn, K. M. 1997, ApJ, 478, 283 
Kenyon, S. J., Calvet, N., & Hartmann, L. 1993, ApJ, 414, 676 



Kurtz, S., Cesaroni, R., Churchwell, E., Hofner, P., & Walmsley, 

C. M. 2000, in Protostars and Planets IV, ed. V. Mannings, A. P. 

Boss, S. S. Russell (Tucson:U. Arizona Press), 299 
Larson, R. B. 2001, in lAU Symp. 200, The Formation of Binary 

Stars, ed. H. Zinnecker & R. D. Mathieu (San Francisco: ASP), 93 
Maxia, C, Tcsti, L., Cesaroni, R., & Walmsley, C. M. 2001, A&A, 

371, 287 

Olmi, L., Cesaroni, R., Hofner, P., Kurtz, S., Churchwell, E., & 

Walmsley C. M. 2003, A&A, 407, 225 
Osorio, M., Lizano, S., & D'Alessio, P. 1999, ApJ, 525, 808 
Osorio, M., D'Alessio, P., MuzeroUe, J., Calvet, N., & Hartmann, L. 

2003, ApJ, 586, 1148 

Rodgers, S. D. & Charnley, S. B. 2001, ApJ, 546, 324 

Rodriguez, L. F., Loinard, L., D'Alessio, P., Wilner, D. J., & Ho, 

P. T. P. 2005, ApJ, 621, L133 
Shepherd, D. S., Claussen, M. J., & Kurtz, S. E. 2001, Science, 292, 

1513 

Scwilo, M., Watson, C, Araya, E., Churchwell, E., Hofner, P., & 

Kurtz, S. 2004, ApJS, 154, 553 
Su, Y.-N., Liu, S.-Y., Lim, J., Ohashi, N., Bcuther, H., Zhang, Q., 

Sollins, P., Hunter, T., Sridharan, T. K., Zhao, J.-H., Ho, P. T. P. 

2004, ApJ, 616, L39 

Terebey S., Shu, F. H., & Cassen, P. 1984, ApJ, 286, 529 
Walmsley, C. M. 1995, Rev. Mexicana Astron. Astrofis. Scr. dc Conf., 
1, 137 

Walsh, A. J., Macdonald, G. H., Alvcy N. D. S., Burton, M. G., & 

Lcc, J.-K. 2003, A&A, 410, 597 
Watson, A. M., de Buizer, J. M., Radomski, J. T., Pifia, R. K., & 

Telesco, C. M. 2003, Revista Mexicana de Astronomia y Astrofisica 

Conference Series, 16, 127 
Watt, S., & Mundy, L. G. 1999, ApJS, 125, 143 
Wilson, T. L., Downcs, D., & Bieging, J. 1979, A&A, 71, 275 
Wood, D. O. S., & Churchwell, E. 1989, ApJS, 69, 831 
Wyrowski, F., Schilkc, P., & Walmsley C. M. 1999, A&A, 341, 882 
Yorkc, H. 1984, Workshop on Star Formation, cd. R. D. Wolstcncroft 

(Edinburgh: Royal Obs.), 63 
Zinchenko, I., Henkel, C, & Mao, R. Q. 2000, A&A, 361, 1079 



Observations and Modeling of MIR Emission from HMPOs 



11 



Table 1 

T-ReCS Filter Information and Integration and Calibration Parameters 





Filter 




Gll.54-0.62 






CM.96-0.02 






G45.07-I-0.13 




Ac 


AA 


t 


Standard 


F„ 


t 


Standard 


F„ 


t 


Standard 


F^ 


2.2 


2.0-2.4 


304 


HD169916 


447 


521 


HD196171 


283 


130 


HD 187642 


550 


3.9 


3.6-4.1 


304 


HD 178345 


67.0 


521 


HD196171 


122 


130 


HD187642 


200 


4.7 


4.4-5.0 


217 


HD 178345 


40.3 


521 


HD196171 


73.5 


130 


HD187642 


137 


7.7 


7.4-8.1 


217 


HD 178345 


17.8 


304 


HD175775 


29.6 


130 


HD169916 


47.7 


8.7 


8.4-9.1 


217 


HD169916 


38.1 


521 


HD169916 


38.1 


130 


HD169916 


38.1 


9.7 


9.2-10.2 


217 


HD169916 


32.3 


304 


HD175775 


20.1 


130 


HD169916 


32.3 


10.4 


9.9-10.9 


217 


HD169916 


28.5 


304 


HD175775 


17.7 


130 


HD169916 


28.5 


11.7 


11.1-12.2 


217 


HD169916 


22.4 


304 


HD169916 


22.4 


130 


HD169916 


22.4 


12.3 


11.7-12.9 


217 


HD169916 


20.0 


130 


HD169916 


20.0 


130 


HD169916 


20.0 


18.3 


17.6-19.1 


217 


HD178345 


3.5 








130 


HD169916 


9.6 


24.6 


23.6-25.5 


304 


Alpha Cen 


27.6 








130 


HD187642 


5.5 



Note: The value 'Af-' is the efreciive central wavelength of the filter in /^m, and 'AA' is the wavelength range of the filter using the 50% transmission 
cut-on and cut-olf of the bandpass. The value 't' is the on-sourcc integration time in seconds in the given filter of the given science held. 'Standard' is 
the name of the standard star used in the given filter for the given science field. 'F,y' is the assumed flux density in Jy of the standard star in the given 
filter. 

" OSCIR data. See De Buizer et al. (2002). 
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Table 2 

Observed flux densities in mJy for sources on Gil. 94-0. 62 field 



A" 


DRT03 1 (HMPO)" 


DRT03 2"= '' 


DR,T03 3"'^ 


DRT03 4"'' 


DRT03 5^'-' 


DRT03 6'^ 


3.9^m 


15±2 


40±5 


71±8 


34±4 


29±3 


52±7 


4.7fim 


37±2 


87±6 


134±8 


106±5 


85±3 


107±9 


7.7/nm 


123±44 


858±308 


2250±808 


1240±446 


854±308 


2420±871 


8.7/jm 


59±5 


505±43 


669±57 


523±45 


433±37 


564±48 


9.7/^m 


14±3 


170±20 


118±17 


120±15 


99±11 


113±18 


10.4;im 


32±4 


412±51 


24()±:i() 


335±42 


296±37 


277±35 


11.7/^m 


94±8 


1210±104 


1010±87 


1150±99 


988±85 


816±71 


12.3jum 


80±19 


1350±316 


1340±316 


1300±305 


1050±247 


1100±258 


18.3/jm 


154±31 


11300±2000 


8800±1570 


11600±2070 


7000±1250 


8230±1470 


2n.8/nn' 


829 ' 340 


")4:in ' 800 


7820 ' 1200 


11300 I IfMO 


11100 ' i5:io 




24.(i;iiii 


296179 


:ir),soui:!20u 


:i:!()00 1 :!010 


-l.'i.lOOl 1800 


2(870012500 





Note: Errors given to the flux values arc 1-a errors that are quadrature addition of the photometric and calibration errors. 

Efi'eetive central filler wavelength. 

Inlcgraled Ilux density in a 0'/6 radius aperture centered on the HMPO candidate. 

This is an extended emission source and is not fully resolved from other sources on the field. 

Integrated flux density in a 2'.'l radius aperture centered on the center of the extended mid-infrared emission from this source. 
^ Integrated flux density in a 2'/ 7 radius aperture centered on the center of the extended mid-infrared emission from this source. 
^ Integrated flux density in a I'.'S radius aperture centered on the center of the extended mid-infrared emission from this source. 
^ Integrated flux density in a l'.'2 radius aperture centered on the center of the extended mid-infrared emission from this source. 

^ This single extended source was labeled as two individual sources in the work of De Buizer et al. (2003, 2005), however the observations presented 
here show it to be a very large and amorphous single but extended source. The integrated flux density is from a 3''l radius aperture centered on the 
center of the extended mid-infrared emission from this source. 
^ From MIRLIN/IRTF data used in De Buizer et al. (2003, 2005). 
Source partially off field at this wavelength so no flux density given. 
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Table 3 

Observed flux densities in mJy for sources on G29.96-0.02 field 



A" 


HMPO" 


HMPO(+bg)'^ 


UC HII" 


2.2/im 


<5 


14 ±2 


327±24 




38 ±5 


112 ±10 


4300±368 


4.7/jm 


123 ±12 


321 ±6 


10800±183 


7.7jum 


66 ±20 


492 ±85 


37400±6440 


8.7Aim 


13 ±2 


437 ±31 


77800±5550 


9.7/jm 


<25 


239 ±26 


56000±6090 


10.4pm 


13 ±3 


370 ±45 


80000±9700 


11.7Mm 


197 ±32 


1130 ±86 


141000±10800 


12.3/im 
18.3/nm'' 


345 ±78 


1680 ±342 


169000±34000 


2280±340 


16400±1640 


462000±46000 



Note: Values preceded hy a means it is a S-a upper limit of the source flux within the apertures specified in the table notes below. Errors given to 
the flux values are i-a errors that are quadrature addition of the photometric and calibration errors. 

Effective central filter vv^avelength. 

Integrated flux density after subtraction of background due to extended flux from the UC H II region. 
^ Integrated flux density in a l''25 radius aperture centered on the HMPO peak without subtraction of the background UC H II region. 

Integrated flux density in a 6''4 radius aperture centered on the center of the extended mid-infrared emission from the UC H II region. 
^ From OSCIR data used in De Buizer et al. (2002). 
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Table 4 

Observed flux densities in mJy for sources on G45.07+0.13 field 



A" 


DRT03 1" 


DRT03 2° 


DRT03 3 (HMPO)" 


DOC05 4' 


DOC05 5f 


DOC05 6" 


DOC05 7" 


2.2/im 


<14 


17±2 


15±2 


<4.04 


g±2 


<4.04 


<4.04 


3.9/xm 


62 ±7 


1460±156 


1440±153 


<3.1 


567±60 


<3.1 


<3.1 


4.7/im 


153 ±10 


4950±81 


4880±80 


<8.3 


1610±27 


<8.3 


<8.3 


7.7/*m 


1330 ±331 


13700±3390 


r750±1910 


72±20 


2710±670 


25±12 


71±20 


8.7([im 


475 ±38 


5940±467 


841±66 


24±3 


996±78 


9±2 


14±2 


9.7Aim 


233 ±27 


2080±197 


165±17 


<14.6 


308±30 


<14.6 


<14.6 


10.4/jm 


377 ±41 


4350±461 


472±50 


17 ±4 


692±73 


9±3 


<9.4 


11.7/jm 


1260 ±88 


18100±1270 


4720±330 


66 ±5 


2350±164 


26±3 


19±3 


12. 3^1111 


2310 ±347 


37100±5560 


11200±1670 


109 ±17 


6610±991 


48±10 


61±11 


18.3/jm 


10900±964 


72900±6360 


12000±1050 


544 ±66 


8600±753 


<136 


264±51 


20.8Mm'' 


15100±2130 


102000'±13500 


19700±2780 










2i.6iJ.Txi 


51700±4150 


564000±45000 


67600±5400 


2170±208 




335±119 


1338±159 



Note: Values preceded by a '<' means it is a 3-a upper limit of the source flux within the apertures specified in the table notes below. Errors given to 
the flux values are l-a errors that are quadrature addition of the photometric and calibration errors. 
Effective central filter wavelength. 

Integrated fiux density in a 2''4 radius aperture centered on center of the extended mid-infrared emission from this source. 

Integrated flux density in a 2'.'0 radius aperture centered on the peak of this source. This source flux density is contaminated by some emission from 
DOC05 5, which is a close mid-infrared companion. DOC05 5 was flt by a gaussian and subtracted off before the aperture photometry was performed 
on DRTOS 2. However DOC05 5 is much fainter than DRT03 2, so the contamination is relatively small (on order the flux error) at all wavelengths. 

Integrated flux density in a I'/l radius aperture centered on the peak of this HMPO candidate. 
^ Integrated fiux density in a 0'/7 radius aperture centered on the peak of this source. 

Flux densities quoted for this source arc rough estimates. DOC05 5 is faint and separated by only 0'/7 from the very prominent mid-infrared source 
DRTOS 2. Fluxes for this source were obtained by fitting DRTOS 2 by a gaussian and subtracting it away before performing aperture photometry with a 
O'.'T radius aperture. However, as the source was observed at longer and longer wavelengths resolution grew steadily worse until at 24.5/Ltm the sources 
could not be resolved. 

^ Integrated flux density in a 0'.'7 radius aperture centered on the peak of this source. 
^ From MIRLIN/IRTF data used in De Buizer et al. (2003, 2005). 

* This is the combined flux from DRTOS 2 and DOC05 5, since they are not resolved from each other using IRTF. 
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Table 5 

Parameters of the Best Fit Models 



HMPO 


r] 






Rc 


L-, 


Pi AU 




i 






(AU) 


(AU) 


(AU) 








n 


Gil. 94-0.62 


2.5 


2 


5000 


30 


75 


1.5x10" 


-13 


53 


G29.96-0.02 


2.5 


245 


12000 


570 


18000 


3.0x10" 


■11 


12 


G45.07+0.13 


2.5 


227 


9000 


370 


25000 


5.3x10" 


12 


35 
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Table 6 

Physical Parameters Derived from the Models 



HMPO 




I'lOOOAu'' 




Spectral'^ 










(Me) 


(K) 




Type 


(Mq yr-i) 


(s-i) 


{Me yr-^) 


Gil. 94-0.62 


1 


40 


5 


AO 


8.7 X 10"^ 


3.0 X 10-^' 


3.5 X 10-^' 


G29.96-0.02 


576 


300 


11 


Bl 


2.0 X 10-2 


1.0 X 10*6 


9.0 X 10-^ 


G45.07+0.13 


64 


200 


13 


Bl 


3.6 X 10-3 


1.0 X 10^'^ 


1.0 X 10-6 



"Mass of the envelope, obtained by integration of the density distribution 
''Temperature of the envelope at a radius of 1000 AU. 

■^Mass of the central star, obtained from L* using the evolutionary tracks of Behrend &c Maeder 2001. 
'^Obtained from tables of Doyon (1990) using the luminosities given in Table 5. 

''Mass accretion rate, obtained from piAU a-nd 

■'^Rate of ionizing photons, obtained from Doyon (1990). 

s Critical mass accretion rate, Merit = 

, where mn is the hydrogen mass, and a is the recombination coefBcient (n 
obtained from Yorke 1984, and Walmsley 1995 assuming spherical accretion. 
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a} b) 




Rfqfit Asceflslon ofTsei (arcsec) Right Ascension offset (ofcsec) 

Fig. 1. — Images of the Gil. 94-0. 62 field taken at a) 4.7 fira, b) 11.7 /im, and c) 18.3 fim at Gemini South with T-ReCS. d) A contour 
plot (smoothed by CK'36) of the 11.7 fim image (blue contours) of the field showing sources with their labels and encircled with dashed lines 
that represent the apertures used for photometry. These aperture radii are given in the notes in Table 2. Plotted as crosses are the locations 
of the water maser groups as given by Hofner & Churchwell (1996). Also plotted (red contours) are the 2 cm radio continuum contours from 
Hofner & Churchwell (1996). The origin of each panel is the mid-infrared peak of the HMPO candidate (R.A.(J2000) = 18'' 14™ 00?29, 
Decl.(J2000) = -18° 53' 23'.' 1). 
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al b) 




10 S D -5 15 10 6 

Right A,5censiori offaet (orcsec), f?lght Aacsnsiori offset (arcaec) 



Fig. 2. — Images of the G29.96-0.02 field taken at a) 4.7 /^m, b) 11.7 /jm, and c) 18.3 fim at Gemini South with T-ReCS. These images had 
to be de-rotated to show north up and east left, hence the gray areas with no data in the corners of the panels, d) A contour plot (smoothed 
by O'.'IS) of the 11.7 fira image (blue contours) of the field showing sources with their labels and encircled with dashed lines that represent 
the apertures used for photometry. These aperture radii are given in the notes in Table 3. Plotted as crosses are the locations of the water 
maser groups as given by Hofner & Churchwell (1996). Also plotted (red contours) are the 2 cm radio continuum contours from Hofner & 
Churchwell (1996). The origin of each panel is the mid-infrared peak of the HMPO candidate (R.A.(J2000) = 18'' 46" 03?74, Decl.(J2000) 
= -02° 39' 22^1). 
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Fig. 3. — Images of the G45.07+0.13 field taken at a) 4.7 ^im, b) 11.7 fim, and c) 18.3 fim at Gemini South with T-ReCS. d) A contour 
plot (smoothed by 0'.'27) of the 11.7 fim image (blue contours) of the field showing sources with their labels and encircled with dashed lines 
that represent the apertures used for photometry. These aperture radii are given in the notes in Table 4. In the upper-right corner there is 
a blow-up of the the region encompassed by the square. This blow-up shows contours of the double source, DRT03 2 and DOC05 5, at 4.7 
fim and is overlaid (red contours) with the 2 cm radio continuum contours from Hofner & Churchwell (1996). Also plotted as crosses are the 
locations of the water maser groups as given by Hofner & Churchwell (1996). The origin of each panel is the mid-infrared peak of the HMPO 
candidate (R.A.(J2000) = 19*" 13™ 22f07, Decl.(J2000) = -1-10° 50' 55'.' 4). 



20 



De Buizer, Osorio, and Calvet 





Fig. 4. — (a) Model SEDs for an = 2 envelope (solid-line) and a TSC envelope (dotted-line) with the same values of the centrifugal radius 
{Rc = 600 AU) and inclination angle {i = 30°). (b) Same as (a), but for a higher value of the inclination angle (i = 60°). (c) Model SEDs 
for two rj = 2 envelopes, with the same value of the inclination angle (i = 30°) and two different values of the centrifugal radius: Rc = 250 
AU (dashed line) and Rc = 600 AU (solid line), (d) Same as (c), but for i = 60°. The stellar luminosity (L* = 25000 Lq), the mass of the 
envelope (Menu = 9 Mq), and the assumed distance (1 kpc) are the same in all models. 
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Fig. 5. — Observed and model SEDs for the source Gil. 94-0. 62 HMPO. The different symbols represent the observed values of the flux 
density. Error bars are S-cr. Upper limits are represented by arrows. Solid-line represents the best fit model, obtained with t = -5.3° , Rc = .30 
AU, L* = 75 Lq, piAU = 2 X 10~^^ g cm""' (sec Table 5), while the dotted line represents a model with the same i and piAU but with 
Rc = 100 AU. Dot-dashed line represents a model with a low inclination angle of i = 30°, with Rc = 30 AU and piAU = 7.5 x lO"^"^ g cm"-^. 
The adopted distance is 4.2 kpc (Hofner & Churchwell 1996). 
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Fig. 6. — Observed and model SEDs for the source G29.96-0.02 HMPO. The different symbols represent the observed values of the flux 
density. Error bars are S-cr for the IR data. The 2.7 mm error bars correspond to the uncertainties assigned by Olmi et al. 2003 due to the 
subtraction method employed. No errors are available for the remaining millimeter flux densities. Upper limits are represented by arrows. 
Solid hne represents the best fit model, obtained with i = f2°. Be = 570 AU, L, = 1.8 X 10* Lq, and piAU = 3 X lO^-'-' g cni^'"* (see Table 
5). Dotted line represents the best fit model obtained by using only the mid-infrared data, with i = 12°, Rc = 380 AU, L* = 2.4 X 10"* Lq, 
and piAU = 5 X 10~^^ g cm~^. Note that the 9.7 /im is only an upper limit and cannot discriminate between the two models in terms of 
the depth of the absorption feature. The adopted distance is 8.4 kpc (Sewilo et al. 2004). Note that both models have the same inclination 
angle, but differ significantly in the luminosity and density. 
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Fig. 7. — Observed and model SEDs for the source G45.07+0.13 HMPO. The different symbols represent the observed values of the flux 
density. Error bars are 3-ct. Upper limits are represented by arrows. Solid line represents the best fit model, obtained with i = 30° , Rc = 370 
AU, Lt = 2.5 X lO** Lq, and pi au = 5-3 x 10^^"^ g cm~^ (see Table 5). Dotted line represents a model with the same luminosity and a higher 
inclination angle, i = 58°, with Rc = 270 AU, and pi au = 7.5 x lO"'^^ g cm"^. The adopted distance is 9.7 kpc (Wood & Churchwell 1989). 



